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ABSTRACT 

The frictional heating by ion-neutral drift is calculated and its effect on 
isobaric thermal instability is carried out. Ambipolar drift heating of one- 
dimensional self-gravitating magnetized molecular slab is used under the assump- 
tions of quasi-magnetohydrostatic and local ionization equilibrium. We see that 
ambipolar drift heating is inversely proportional to density and its value in some 
regions of the slab can be significantly larger than the average heating rates of 
cosmic rays and turbulent motions. The results show that the isobaric thermal 
instability can occur in some regions of the slab; therefore it may produce the 
slab fragmentation and formation of the AU-scale condensations. 

Subject headings: ISM: clouds - ISM: structure - instabilities - star: formation 

1. Introduction 

Contour maps of the molecular clouds show a very fragmented distribution, consisting 
several clumps and cores (see, e.g., Williams et al. 2000, Falgarone et al. 2004). Cloud 
structures are likely to be more complex than this because many of the clumps that are 
optically thick in one transition contain substructure that is visible in other transitions 
(Perault et al. 1985, Kwan & Sanders 1986, Dickman et al. 1990). 

The reality of AU-scale structure within molecular gas is still questionable. Indeed, only 
minor species can be easily observed towards molecular clouds and their spatial distribution 
might not reflect that of H 2 were present. H 2 CO, HCO + , and OH apparently display column 
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density fluctuations reaching 5 to 15% along lines of sight separated by 10 AU (Marscher 
et al. 1993, Moore & Marscher 1995, Liszt & Lucas 2000), while dust grains appear to be 
more smoothly distributed (Thoraval et al. 1996). Rollinde et al. (2003) probed the spatial 
distribution of CH, CH + , and CN at scales ranging from 1 AU up to 150 Au. Studies of the 
time variability of absorption lines indicates the presence of fluctuations on scales of 10~ 4 pc 
(5-50 AU) and masses of 1O~ 9 M (Boisse et al. 2005). At larger scales (about 10,000 AU), 
Pan et al. (2001) find significant differences in CN, CH, and CH + absorption lines. 

Star forming regions are often interpreted in a picture that consists of an ensemble 
of discrete regions of gas. Existence of these AU-scale condensations in the clumps of the 
molecular clouds, may be a motivation to investigate the hierarchy down-up structure of 
the clouds. Various authors, such as Kwan (1997) and Silk & Takahashi (1997), have cited 
coagulation processes as an origin for the form of the initial mass function (IMF). For exam- 
ple, Murray & Lin (1996) have argued that in protoglobular clusters, gas is fragmented onto 
sub- Jeans-mass clumps by a variety of hydrodynamical instabilities, and that eventual star 
formation ensues when enough clumps have merged to form Jeans unstable objects. Pringle 
et al. (2001) have considered a similar process in the assembly of giant molecular clouds in 
the spiral arms of disc galaxies. Thus, the condensations appear, as a possible scenario, to 
be precursors of large-scale cores (dense cores with significant jeans mass) via merging and 
collisions; they constitute the initial conditions for star formation. This idea is largely an 
attempt to characterize the observed structure of molecular cloud complexes- see, e.g. the 
review by Williams et al. (2000)- but is likely to be applicable in many situations. 

It is well known that molecular clouds are turbulent which triggers the formation of 
structures. A recent literature has suggested that the dumpiness in clouds arises naturally 
from their formation through thermal instability, which acts on time-scales that can be 
much shorter than the duration of the turbulent motions. This suggests that molecular 
clouds may be already fragmented when they form (Koyama & Inutsuka 2002, Audit & 
Hennebelle 2005, Heitsch et al. 2006, Vazquez-Semadeni et al. 2006). In the same way, other 
mechanisms to generate substructures inside molecular clouds, based on MHD waves, have 
been proposed (Gammie & Ostriker 1996, Falle & Hartquist 2002, Folini et al. 2004). In 
this paper we make the assumption that the molecular cloud is initially an uniform ensemble 
which then fragments due to thermal instability. This is clearly a competitive mechanism 
that we considered although all could operate. 

It is also believed that the interstellar medium is strongly magnetized. Although the 
structure of the magnetic field is poorly known, various studies have considered thermal 
instability and magnetic field as well. Loewenstein (1990) investigated the effect of magnetic 
fields on thermal instability in cooling flows using linear perturbation analysis. Steele & 



- 3- 



Ibahez (1999) consider the nonlinear thermal disturbances for a two-dimensional structure 
taking into account thermal conduction parallel to and perpendicular to the magnetic field. 
Hennebelle & Perault (2000) have investigated the role of an initially uniform magnetic 
field, analytically and numerically, when the thermal conduction is dynamically triggered. 
They propose a mechanism based on magnetic tension to explain the thermal collapse in 
a magnetized flow and argue that the magnetic intensity in the WNM and in the CNM 
should not be very different. The effect of Alfven waves on thermal instability was studied 
by Hennebelle & Passot (2006). Inoue et al. (2007) Performed one-dimensional two-fluid 
numerical simulations that describe the neutral and ionized components in the interstellar 
medium with purely transverse magnetic fields. They find that the cloud that is formed 
by the thermal instability always has magnetic field strength on the order of a few fiG 
irrespective of the initial strength. 

In molecular clouds, a neutral molecular gas intermixed with an ionized component 
that is tied directly to the magnetic field. Birk (2000) has used the two-fluid technique to 
find the thermal condensation modes in weakly ionized hydrogen plasma. Nejad-Asghar & 
Ghanbari (2003) studied the effect of linear thermal instability in a weakly ionized magnetic 
molecular cloud within the one-fluid description. They concluded that not only the thermal 
instability is an important mechanism but also there are solutions where the isentropic ther- 
mal instability including ambipolar diffusion allows compression along the magnetic field but 
not perpendicular to it. Stiele et al. (2006 hearafter SLH) revisit carefully the above sce- 
nario with some improvements. They drive a critical wavelength that separates the ranges of 
stability and instability. Moreover, they discussed on the relations between the parameters 
of cooling and heating rates, where isobaric and/or isentropic thermal instability can be oc- 
curred. Nejad-Asghar & Ghanbari (2006) investigated the effect of self-gravity in occurrence 
of the thermal instability in a weakly ionized magnetic molecular cloud. They concluded 
that there are, not only oblate condensation forming solutions, but also prolate solutions 
according to expansion or contraction of the clump. 

The goal of this paper is to analyze the effect of ambipolar drift heating on net cooling 
function and investigation of the isobaric thermal instability in a self-gravitating magnetized 
molecular slab. We re-formulate the equations of a quasi-magnetohydrostatic self-gravitation 
slab in § 2. Then, density dependence of the drift velocity is investigated. Section 3 deals 
with the net cooling rate of a standard molecular cloud as discussed by Goldsmith (2001, 
hereafter G01). Moreover, we investigate occurrence of the isobaric thermal instability and 
formation of the AU-scale condensations in some regions of the molecular slab. Finally, § 4 
allocates to a gathering summary with conclusion. 
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2. Self-Gravitating Slab 

We consider the evolution of a plane-parallel self-gravitating magnetized slab of lightly 
ionized molecular gas (see, Fig. [1]). In this geometry, all variables are functions of distance 
z to the central plane and time t only. The magnetic field is frozen only to the ions, so that 
diffusion of the field relative to the neutral gas must continuously insert a drag force (per 
unit volume) 

fd = lADpipVd, (1) 

where 'Jad ~ 3.5 x 10 10 m 3 /kg.s represents the collision drag coefficient and pi = ep u is the 
ion density. In ionization equilibrium state v = 1/2 and e ~ 9.5 x 10 -15 kg 1//2 .m -3 / 2 is valid, 
and v ^ 1/2 represents a deviation from equilibrium state (Shu 1992). The drift velocity is 
given by 

1 d . B\ 

which is obtained by assumption that the pressure and gravitational force on the charged 
fluid component are negligible compared to the Lorentz force because of the low ionization 
fraction. In this way, the magnetic fields are directly evolved by charged fluid component, 
as follows: 

dB dv d , 

Since the ion density is negligible in comparison to the neutral density, the mass con- 
servation is given by 

dp dv 

dt = - P d~z- (4) 

The momentum equation then becomes 

^ V — g 1 ^ (a 2 p + ^ ) (5) 
dt p dz 2p 

where a is the isothermal sound speed and the gravitational acceleration g obeys the poisson's 
equation 

Yz = - AnGp - (6) 

Following the many previous treatments, a further simplification is possible if we intro- 
duce the surface density between mid-plane and z > as 



a= / p(z',t)dz'. (7) 
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The usual denned surface density from — z to +z has twice the value of a. By transformation 
from (z, t) to (a, t), 



d ( d \ d f d 



dt \dt I dz ^ V da 



the drift velocity is 



1 d B 2 , 
^ADep y da 2/i ' 
and the equation becomes 



d ,B, Id ,B 2 dB , 



(9) 



dt p 1ad£Po da p v da 
With the above, field equation ([6]) can be integrated to give 

g = -AnGa, (10) 

while the equation of continuity (HI) and the equation of motion ((5]) take the form 

dz _ 1 
da p 

and 

d 2 z . „ /2 B 2 



11) 



respectively. The slab is assumed to be in quasi-magnetohydrostatic equilibrium at all time, 
supported against its own self-gravity by the magnetic nd gas pressure. The loss of flux 
from ambipolar diffusion is exactly compensated for by the compression of the slab which 
is necessary to maintain equilibrium. In this approximation, the left-hand side of equation 
(|T2|) is zero, and we may integrate the force balance to obtain 

f- + a 2 p = 2rrG(al - a 2 ) (13) 

where integration constant a^ is the value of a at z = oo (where p is zero). 

In a sense, our goal is to obtain the variation of drift velocity as a function of density 
at initial time. Evolution of the slab is beyond the scope of this paper. For this purpose, we 
follow the work of Shu (1983) to introduce the non-dimension quantities 



a 



2 _ 2nGa 



oo 



° = — ' P = P» Z = 15 Z i 



(Too ' 2nGal 



a 



b s -^, N^K^K d4) 
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so that we rewrite the basic equations (jHJ), (fTTl) and ([TBI as follows: 

9 ( B d B 2 dB 



(15) 



£ 2 + p = l-a 2 , (16) 

P = i (17) 

In this way, a natural family of initial states is generated by assuming that the initial ratio 

of magnetic to gas pressure is everywhere a constant, a>o, i.e., B 2 /p = «o at i = 0. Then one 
finds from equations (|T6)) . (ITT)) and (|HD that 

Po /1oN 

P(z,t=o) 



cosh (z/zoo) 

„, = „ 2a ° f^y\AT^>, (i 9 ) 



VI + «o 7ADePo V2 V Po / V Po 
where p = 27rG<7 2 XD /a 2 (l + a ) is the central density of the slab at t — 0, and z^ = 
a-^/ (1 + a )/2nGp is a length-scale parameter. 

The interior magnetic field of a molecular cloud acts only on the ions and electrons, so 
that the tendency for the field lines to straighten out is impeded by collisions of charged 
particles with the neutral particles. A simple physical derivation of the ambipolar drift 
heating rate (J.kg _1 .s _1 ) is given by considering the drag force (ED) as follows 

Tad = — = l AD ep v vl (20) 

Substituting the equations ([TBI) and (jT9)) into equation ([20)) . the ambipolar drift heating rate 
of a molecular slab is given by 

qg 8nG 2 ^ 1 - p/p 

i -AD = — a Po -f—, — (21) 

l + «o7ADe (P/Po) 

which can be compared to the heating of cosmic rays (e.g., Goldsmith 2001), 

T CR « 3.12 x 10" 8 J.kg _1 .s _1 , (22) 

as follows (y — 1/2) 

Tad _ ®l { a \ 2 ( n o \ x l 2 \ — n/nQ 



T C r l + a VSOOm.s" 1 / Vl0 12 m- 3 / {n/n Q ) l l 
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that is shown in Fig. [2] for typical values of «o = 1, a, = Cl.3km.s~ 1 , and no = 10 12 m -3 . 
We have chosen to plot the density profile ([TBI , because this makes it easy to compare the 
relative magnitude of ambipolar drift heating at different regions of the slab. Since the drift 
velocity is inversely proportional to the density (e.g., Eq IT9~j) . it is clear that the ambipolar 
drift heating rate is predominated in the outer regions of the slab. 

It would also be interesting to derive a rough estimate of heating due to the drift induced 
by turbulent motions and to compare its order of magnitude with the heating of cosmic rays. 
Following Black (1987) the turbulence dissipation heating rate can be estimated as 

r {l mv turb)( v turbn)(j) 
1 turb ~ 

mn 

- 1.6 x io- 8 (-^r) 3 (i^) J.kg-i.s- 1 , (24) 
1 km.s t 

where v tur b is the turbulent velocity and / is the dissipation length scale. With v tur b ~ 
1 km.s -1 and I ~ 1 pc, we calculate T tur i, ~ 1.6 x 10 -8 J.kg -1 .s -1 , comparable to half of the 
heating rate of cosmic rays (i.e., Eq. I2"2"j) . 



3. Isobaric Thermal Instability 

3.1. The Molecular Gas Cooling Function 

The cooling rate for molecular gas that is used in this work is based on the "standard 
abundances" adopted by G01. I will briefly review here the important aspects of its construc- 
tion. Calculation of thermal balance in molecular clouds by Goldsmith & Langer (1978), 
Neufeld et al. (1995), and G01 have shown that the various isotopic variants of the carbon 
mooxide molecule are certainly important coolants at high densities (i.e. n > 10 10 m -3 ). 
This is because the transitions that contribute to the cooling are thermalized due to the 
combination of higher collision rate and radiative trapping. The "standard abundances" 
for 12 CO, 13 CO, and C 18 was adopted by G01 based on the observational data equal to 
5.62 x 10" 5 , 1.00 x 10~ 6 , and 1.00 x 10~ 7 relative to that of H 2 , respectively. Neutral Carbon 
is an important coolant since although it has only two fine structure transitions available, 
they have modest energies, and the abundance of this atom is relatively large. 

The fractional abundances of ortho and para water vapor in molecular clouds are sur- 
prisingly low, so that G01 adopted a number density of ortho-H 2 relative to that of H 2 
equal to 10~ 8 . He have also doubled the ortho-H 2 cooling rate to account approximately for 
contribution of para-H 2 0. There is a large number of molecular species, including CN, HCN, 
CS, NH 3 , HCO + , and H 2 CO, that have low fractional abundances. To simplify the cooling 
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calculation, G01 proposed the CS molecule as representative and multiplied its cooling by a 
factor of 10 to include the contribution of the other species in this category. 

The velocity field in molecular clouds is not well understood in detail. We can assume 
that the large velocity gradient approximation is valid for the radiative transfer calculations 
necessary to compute the molecular cooling when the optical depths are not small. A critical 
input to the large velocity gradient calculation is the molecular abundance per unit velocity 
gradient, which is equivalent to the column density per unit line width. Comparing the sizes 
and line widths for many clouds gives a characteristic value of 1 km.s _1 .pc _1 that was used 
by G01. The use of the large velocity gradient model introduces some uncertainty, but, do 
not significantly alter the total cooling rate. 

In this way, the total cooling has a complicated dependence on the temperature and 
density of molecular hydrogen, parameterized by G01 as 

A(n,T) = A(n)(^) /3{n) J.S^.kg- 1 , (25) 

where T = 10K and the values of A(n) and (3(n) are given in Table 1 for standard abundances 
and velocity gradient of 1 km.s -1 .pc -1 . Here, we use the polynomial fitting functions as 
follows 

l0g ( iW-iVi ) = " °- 87(1 ° g V ] ~ °- 14(1 ° g ^ )2 ' (26) 

\ J. Kg .S J Uq Uq 

(3 {n) = 3.07 - O.llQog -) - 0.13(log -) 2 , (27) 
where no = 10 12 m~ 3 . Another expression for the parameter A( n ) is 

Tl 

logA (n) =logA + 51og— , (28) 

n 

which was used by Nejad-Asghar & Ghanbari (2003, 2006) and SLH, for linear investigation 
of thermal instability in a narrow range of density. In this formalism, the approximate values 
of A and 5 between neighbor densities are given in Table 1. 

There is observational evidences that a variety of molecular species are depleted on the 
surface of dust grains (e.g., Tielen & Whittet 1997). The observational results on depletion of 
different species are relatively limited and are certainly compromised by blending of emission 
from different regions along the line of sight. It is quite difficult to know exactly how the 
column densities or abundances of species observed on dust grain surfaces are related to the 
abundances in the gas phase. There also have not been systematic studies of the depletion of 
a wide range of molecular species in a set of sources. In this way, G01 used the fragmentary 
data from a selection of clouds to analyze the effects of molecular depletion on the standard 
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values of the cooling function parameters A(n) and /3r n )- Because of the above mentioned 
uncertainties, we should treat the cooling function in a parametric fashion with A(n) and 
/3( n ) approximately near the standard values of Table 1. 

The cooling function was constructed assuming that statistical steady-state equilibrium 
is valid, so that the chemical times are very short with respect to dynamical and cooling times. 
For all of the trace coolants this is reasonable except H 2 molecule that the collision time 
(time to reach statistical steady-state equilibrium) is necessarily of order the cooling time 
(Gilden 1984). The H 2 is the main coolant at higher temperatures (the temperature where 
the carbon monoxide and H2 contribute equally is slowly increasing function of density). 
At the very low temperatures of interest considered by G01 in the molecular clouds, H2 
do not play a significant role as a major coolant because of its low-lying energy levels. 
Thus, the standard values of the cooling function parameters A(n) and f3t n ) (as mentioned 
in Table 1) based on the statistical steady-state equilibrium are approximately appropriate 
for our analysis. 

Models of the molecular clouds identify several different heating mechanisms. Here, 
we use the heating of cosmic ray (1221) and the ion-neutral slip heating rate as described by 
equation (1231) . The isobaric instability mode in a thermal equilibrium gas can be discussed 
in relation to the two-phase model in the same way as described in figure 3 of Gilden (1984). 
The detailed study of the linear investigation and dynamical evolution of the isobaric modes 
are given in the next sub-section. The instance here can be analyzed in a graphical approach. 
The pressure-density plane is an elegant method to describe the isobaric instability modes 
in the thermal equilibrium, 

A {n) (^r {n) =r C R + r AD , (29) 
-to 

of a molecular gas. Substituting the equation ( 1231 into ( 1291) . the pressure in thermal equi- 
librium self-gravitating molecular slab can be obtained as 



p nT 



Po n T 

n , f Y CR 



n { A(n) 



1 + 



2 1/11 l//3(«) 



1 + ct 300m.s- 17 10 12 m~ 3/ (n/n ) 1/2 . 



(30) 



The two-phase instability is displayed in Fig. [3] where the values of a and no are typically 
chosen equal to 0.3km.s _1 and 10 12 m -3 , respectively. The diagram of Fig. [3] shows the 
equilibrium state without magnetic field (i.e., «o = 0) in comparison with considering of 
the ambipolar drift heating rate with «o = 10. This figure suggests that without ambipolar 
diffusion friction, there would not be isobaric thermal instability. This is not what has been 
found by Gilden (1984), because as we known the thermal instability is very sensitive to the 
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net cooling function. Here we use the approximately most accurate cooling function which is 
evaluated by GO. An equilibrium state is specified by the intersection of a thermal equilibrium 
curve and a constant pressure line. The two curves in Fig. [3] display that considering of 
ambipolar drift heating rate causes to bifurcate the slab into two phases, with individual gas 
elements going into that phase corresponding to the sign of the initial fluctuation. 



3.2. Condensation Growth 

Equation ( |T9l) shows that the drift velocity is inversely proportional to the density so 
that we can choose, in a general form, = const, x p~ b where b > 0. For a self-gravitating 
slab, the value of b may be approximated roughly between 0.5 to 2.0. This range of values 
should be compared to equation (10) of Nejad-Asghar & Ghanbari (2006), where b = 1 + v 
was used. The heating due to ion-neutral friction is spatially dependent (because it depends 
on gas density and gas density varies with position). Here we do not consider this issue 
in the instability analysis because we are interested considering the AU-scale condensations 
that may be formed from waves with small wavelengths (respect to the cloud size). In this 
formalism, equation (37) of SLH for isobaric thermal instability may be generalized to 

A cl = , A ° (31) 

with £ = r ad/ 'A and Ao = 2n^KT/ p/3A where K is the thermal conduction coefficient. If 
the wavelength of perturbation is larger than the critical wavelength (13 ip . there will occur 
isobaric instability. This requires 

5 + (2b - v)£ > 13. (32) 

Since < £ < 1, we find that 

5> (3-(2b-u) (33) 

is required for occurrence of isobaric thermal instability. On the other hand, critical wave- 
length of isentropic thermal instability (equation 38 of SLH) is modified as 

A c2 = , A ° (34) 

/ 3 8-(2b-u)j 7 
V"2 fi 1 

which requires 

5<~p-(2b-u)t (35) 
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Thus, the condition 

5 < (36) 

must at least be held for occurrence of isentropic thermal instability. The instability regimes 
are summarized in Fig. HI Depending on the chosen density range, the values of (3 and 5 can 
differ, but we can only expect the isobaric thermal instability mode. The isentropic regime 
can not be fulfilled in the molecular slab. Therefore, we only consider the nonlinear isobaric 
thermal instability mode. 

To give a complete description of the nonlinear thermal runaway, one needs to solve 
the nonlinear fluid equations. The energy equation follows from the use of the first law of 
thermodynamics, that is 

1 dp 7 Pdp , n n 

—Tt-—r P Tt +p n <"*> = ' (37) 

where 7 is the polytropic index, fi( Pi T) is the net cooling function 

^(p,t) = A( Pi t) - (T C r + Tad), (38) 
and the pressure p is given by the ideal gas equation of state 

V = ~pT, (39) 
P 

where fi is the mean atomic mass per particle and R is the universal gas constant. Being 
interested in the development of the isobaric thermal instability in the quasi-hydrostatic 
molecular slab, we put dp/dt = in the energy equation (1371) and obtain the following 
equation 

where the isobaric condition p oc T _1 is used. Now let us introduce the dimensionless 
temperature, T = T/Tq, to present the energy equation (HOI as dimensionless form 



dT 

~dt 



Y + «(f A = (41) 



where Ct/f s \ is the dimensionless net cooling function of the slab 

Q ™ = { ^RTj { 2(2nGy^ ){ 2 7 rGa^ )n ^- (42) 

Thus, the nonlinear dynamics of the isobaric thermal instability is determined solely by the 
form of the net cooling function curve Q,f « and the initial conditions in the equation 
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We use the equations (j2Tj) . fl22|) . and fl25l) - fl27j) to find the net cooling function of a self- 
gravitating molecular slab. Figure [5] shows the net cooling function of the slab for typical 
values of T^ =0 y The necessary condition for the instability has the form dCl/dn > 0, which 
in the slab can be expressed as dfl/dz < 0. This condition can be easily seen to be fulfilled 
with the equilibrium points u\ and v>2 in the Fig. [5j Indeed, these results accurately coincide 
with the linear instability regimes that are summarized in Fig. [3j We consider a constant 
initial form of the temperature profile in the slab, and solve numerically the temperature 
equation (141 p with approximation that the dynamics of the slab is solely determined by the 
equations f|T5|) - f|T7|) and the initial condition ffTB"]) as outlined by Shu (1983). Figure [6] gives 
a plot of the evolution of the temperature profile in a typical case described in Fig. [51 As 
can be seen in Fig. [61 large temperature gradients are created with time at the interfaces 
between different unstable thermal phases. 



4. Summary and Conclusion 

In this paper the importance of ambipolar drift heating rate on net cooling function 
have been investigated, and the isobaric thermal instability in a self-gravitating magnetized 
molecular slab is analyzed. We re-formulate the equations of a quasi-magnetohydrostatic 
self-gravitation slab so that the density dependence of the drift velocity is estimated. Since 
the drift velocity is inversely proportional to the density, the ambipolar drift heating rate is 
predominated in the outer regions of the slab. 

Choosing a parametric density dependence of drift velocity, the isentropic and isobaric 
thermal instability criteria of SLH are generalized. We have applied the cooling rate of G01 
to find the physical contents of the instability criteria. We have found that the isentropic 
thermal instability can not occur in the molecular cloud, while the isobaric criteria is ac- 
curately complied. This term is also investigated by considering of thermal equilibrium in 
pressure-density plane, and by plotting of the net cooling function in Fig. [51 

We have considered a constant initial form of the temperature profile in the slab, and 
have solved numerically the nonlinear isobaric energy equation with approximation that 
the dynamics of the slab is solely determined by the isothermal equations as outlined by 
Shu (1983). In this way, a large temperature gradient are created with time at the interfaces 
between different isobaric unstable thermal phases. In the context of our thermal instability 
analysis in the molecular slab, we find that the isobaric thermal instability can occur in some 
regions of it, therefore it may produce the slab fragmentation and formation of the spherical 
AU-scale condensations. 
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Fig. 1. — Schematic diagram of self-gravitating magnetic slab in z direction. 
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Fig. 2. — Relative magnitude of the ambipolar drift heating rate of a self-gravitating molec- 
ular slab in the case of a,® = 1, a — 0.3km.s _1 , and n = 10 12 m~ 3 . The density profile of the 
slab is given for comparison. 



-17- 




Fig. 3. — Curves of thermal equilibrium in the pressure-density plane. The dash curve 
represents the equilibrium state without magnetic field (i.e., a = 0), and the solid curve is 
obtained by considering the ambipolar drift heating rate with «o = 10, a = 0.3km.s _1 , and 
n = 10 12 m~ 3 . 
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Fig. 4. — Isobaric (horizontal shade) and isentropic (vertical shade) thermal instability 
regimes in terms of temperature dependence (3 and density dependence 6 in the parameterized 
cooling function. The squares □ represent the values of (3 and 5 from table 1. 
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Fig. 5. — The non-dimensional net cooling function of the self-gravitating slab with ao = 10, 
a = 0.3km.s _1 , and n = 10 12 m~ 3 for cases Ta =0 \ = 4 (dash), 5 (solid), and 6 (dot). The 
points «i and u 2 are unstable because of dQ/dz < 0. 



-20- 




0.0 0.5 1.0 1.5 2.0 



Fig. 6. — Evolution of the temperature of the isobaric unstable regions in self-gravitating 
molecular slab with «o = 10, a = 0.3km.s _1 , and n = 10 12 m~ 3 from initial value T = 5 at 
t = (dash) to a typical time t = 1 (solid). 
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Table 1: Parameters 3, for the gas cooling function with standard abundances 13 and velocity 
gradient of 1 km.s~ 1 .pc~ 1 . 



-4 


-7.82 


1.4 


}^ 


-6.14 


0.42 


-3.5 


-7.61 
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}^ 


-7.40 


0.06 


-3 


-7.58 


2.4 


}^ 


-8.45 


-0.29 


-2 


-7.87 


2.7 


}^ 


-8.65 


-0.39 


-1 


-8.26 


3.0 


}^ 


-8.96 


-0.7 





-8.96 


3.4 









a A (n , T) =A ( „ ) (T/10K)' 3 <"). 
^Goldsmith (2001). 
c n = 10 12 irT 3 . 
d A (n) = Ao(n/n ) 5 . 



